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1. Executive Summary 

There is accumulating evidence for a difference between neutrino and antineu- 
trino oscillations at the ~ 1 eV^ scale. The MiniBooNE experiment observes an 
unexplained excess of electron-like events at low energies in neutrino mode^^, which 
may be due, for example, to either a neutral current radiative interaction^^^terile 
neutrino decay or to neutrino oscillations involving sterile neutrinos SEEITIHini 

rfni 

and which may be related to the LSND signal^. No excess of electron-like events 
0.5 ± 7.8 ± 8.7), however, is observed so far at low energies in antineutrino mode 
. Furthermore, global 3-1-1 and 3-1-2 sterile neutrino fits to the world neutrino and 
antineutrino datal^ suggest a difference between neutrinos and antineutrinos with 
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significant (sin^ 26^^^ ~ 35%) disappearance. In order to test whether the low- 
energy excess is due to neutrino oscillations and whether there is a difference between 

and disappearance, we propose building a second MiniBooNE detector at (or 
moving the existing MiniBooNE detector to) a distance of ~ 200 m from the Booster 
Neutrino Beam (BNB) production target. With identical detectors at different dis- 
tances, most of the systematic errors will cancel when taking a ratio of events in the 
two detectors, as the neutrino flux varies as to a calculable approximation. This 
will allow sensitive tests of oscillations for both z/g and Ve appearance and and 

disappearance. Furthermore, a comparison between oscillations in neutrino mode 
and antineutrino mode will allow a sensitive search for CP and CPT violation in the 
lepton sector at short baseline (Am^ > 0.1 eV^). Finally, by comparing the rates for 
a neutral current (NC) reaction, such as NC 7r° scattering or NC elastic scattering, a 
direct search for sterile neutrinos will be made. The initial amount of running time 
requested for the near detector will be a total of ~ 2E20 POT divided between neu- 
trino mode and antineutrino mode, which will provide statistics comparable to what 
has already been collected in the far detector. A thorough understanding of this 
short-baseline physics will be of great importance to future long-baseline oscillation 
experiments. 



2. Introduction 

Evidence for neutrino oscillations comes from solar-neutrino I 13 | 14 | 15 | 16 | 17t 
reactor-antineutrino experiments ^1^, which have observed z/g disappearance at Am^ ~ 
8 X 10~^ eV^, and a tmosp heric-neutrino 19 20 2122) long-baseline accelerator- 
neutrino experiments I23l22i^ which have observed disappearance at Am^ ~ 3xl0~^ 
eV^. In addition, the LSND experiment has presented evidence for z/^ — > z/g os- 
cillations at the Am^ ~ 1 eV^ scale. If all three phenomena are caused by neutrino 
oscillations, these three Am? scales cannot be accommodated within an extension of 
the Standard Model with only three neutrino mass eigenstates. An explanation of 
all three mass s cales with neutrino oscillations requires the addition of one or more 
sterile neutrinos I 4 | 5 | 6 | 7 |819ji further extensions of the Standard Model {e.g.,^^). 

The MiniBooNE experiment was designed to test the neutrino oscillation inter- 
pretation of the LSND signal in both neutrino and antineutrino modes. MiniBooNE 
has approximately the same L/Ey as LSND but with an order of magnitude higher 
baseline and energy. Due to the higher energy and dissimilar event signature, Mini- 
BooNE systematic errors are completely different from LSND errors. MiniBooNE's 
updated oscillation results in neutrino mode ^ show no significant excess of events 
at higher energies; however, a sizeable excess of events is observed at lower ener- 
gies. Although the excess energy shape does not fit simple two-neutrino oscillations, 
the number of excess events agrees approximately with the LSND expectation. At 
present, with 3.4E20 POT in antineutrino mode, MiniBooNE observes no excess so 
far at lower energies, while at higher energies the data are inconclusive with respect 
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to antineutrino oscillations suggested by the LSND data'^. 

A global fit to the world neutrino oscillation data has been performed in terms of 
3+1 and 3+2 sterile neutrino models'^. Although a large incompatibility is observed 
when fitting all of the data, there is much less tension when fitting the world neutrino 
and antineutrino data separately. Figs. [H El [3l and H] show the allowed regions 
from global fits to the world neutrino and antineutrino data, assuming 3+1 sterile 
neutrino models. The best neutrino fit occurs at Amli = 0.19 eV^ with sin^ 26*^^6 = 
0.031, sin^2^^^ = 0.031, and sm'^2eee = 0.034, corresponding to a = 90.5/90 
DF (probability = 47%). The best antineutrino fit occurs at Amli = 0.915 eV^ 
with sin^ 26*^2 = 0.0043, sin^ 26*^^ = 0.35, and sin^ 26'ee = 0.043, corresponding to a 

= 87.9/103 DF (probability = 86%). The antineutrino fit is dominated by the 
LSND data^^; however. Fig. [3] also shows that the global antineutrino data without 
LSND is consistent with the LSND signal and has a closed 90% CL contour around 
the LSND best-fit point. CPT-conserving oscillation scenarios appear insufficient to 
explain all of the data. As stated in reference "CPT- violating oscillations or 
effective CPT violation I^HM may succeed in reconciling all short-baseline oscillation 
signatures, and should be pursued". 

3. MiniBooNE 



3.1. Description of the Experiment 

A schematic drawing of the MiniBooNE experiment at FNAL is shown in Fig. O 
The experiment is fed by 8-GeV kinetic energy protons from the Booster that interact 
in a 71-cm long Be target located at the upstream end of a magnetic focusing horn. 
The horn pulses with a current of 174 kA and, depending on the polarity, either 
focuses TT"*" and and defocuses vr" and K~ to form a pure neutrino beam or 
focuses 7r~ and K~ and defocuses 7r+ and to form a somewhat pure antineutrino 
beam. The produced pions and kaons decay in a 50-m long pipe, and a fraction of the 
neutrinos and antineutrinos interact in the MiniBooNE detector, which is located 
541 m downstream of the Be target. For the MiniBooNE results presented here, a 
total of 6.5 X 10^° POT were collected in neutrino mode and 3.4 x 10^° POT were 
collected in antineutrino mode. 

The MiniBooNE detector consists of a 12.2-m diameter spherical tank filled 
with approximately 800 tons of mineral oil {CH2). A schematic drawing of the Mini- 
BooNE detector is shown in Fig. [61 There are a total of 1280 8-inch detector pho- 
totubes (covering 10% of the surface area) and 240 veto phototubes. The fiducial 
volume has a 5-m radius and corresponds to approximately 450 tons. Only ~ 2% of 
the phototube channels failed over the course of the run. 
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Figure 1: The allowed region for appearance from a global fit to the world neutrino data, assum- 
ing a 3+1 sterile neutrino model. The star indicates the best-fit point at Amli = 0.19 eV^ and 
sin^ 20ne = 0.031. 
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Figure 2: The allowed region for disappearance from a global fit to the world neutrino data, 
assuming a 3+1 sterile neutrino model. The star indicates the best-fit point at Am^i — 0.19 eV^ 
and sin^ 26 nn = 0.031. 
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Figure 3: The allowed region for appearance from a global fit to the world antineutrino data, 
assuming a 3+1 sterile antineutrino model. The star indicates the best-fit point at Arn^i = 0.915 
eV^ and sin^ 20^e = 0.0043. Also shown is the global fit to the world antineutrino data without 
LSND. 
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Figure 4: The allowed region for disappearance from a global fit to the world antineutrino data, 
assuming a 3+1 sterile antineutrino model. The star indicates the best-fit point at Am^^i = 0.915 
eV^ andsm'^20nu =0.35. 



7 




MiniBooNE Detector 



Signal Region 




Figure 6: A schematic drawing of the MiniBooNE detector. 
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Figure 7: The v^^ CCQE distribution for data (points with error bars) compared to the MC 
simulation (histograms). 



3.2. MiniBooNE Cross Section Results 

MiniBooNE has published two cross section results. First, MiniBooNE has made 
a high-statistics measurement of charged- current quasi-elastic (CCQE) scattering 
events Fig. [7] shows the CCQE distribution for data (points with error 
bars) compared to a MC simulation (histograms). A strong disagreement between the 
data and the original simulation (dashed histogram) was first observed. However, by 
increasing the axial mass, M^-, to 1.23 ±0.20 GeV and by introducing a new variable, 
K = 1.019 ± 0.011, where k is the increase in the incident proton threshold, the 
agreement between data and the simulation (solid histogram) is greatly improved. It 
is impressive that such good agreement is obtained by adjusting these two variables. 

MiniBooNE has also collected the world's largest sample of neutral-current 7r° 
events as shown in Fig. [81 By fitting the 77 mass and Et^II — cosOt^) distributions, 
the fraction of vr^ produced coherently is determined to be 19.5±1.1±2.5%. Excellent 
agreement is obtained between data and MC simulation. 

3.3. Neutrino Oscillation Event Selection 

MiniBooNE searches for —>■ z/g oscillations by measuring the rate of u^C —>■ e~X 
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Figure 8: The neutral- current tt*^ 77 mass and i?Tr(l — cos^tt) distributions for data (points with 
error bars) compared to the MC simulation (histograms). 



CCQE events and testing whether the measured rate is consistent with the estimated 
background rate. To select candidate u^. CCQE events, an initial selection is first 
applied: > 200 tank hits, < 6 veto hits, reconstructed time within the neutrino 
beam spill, reconstructed vertex radius < 500 cm, and visible energy E^s > 140 
MeV. It is then required that the event vertex reconstructed assuming an outgoing 
electron and the track endpoint reconstructed assuming an outgoing muon occur at 
radii < 500 cm and < 488 cm, respectively, to ensure good event reconstruction 
and efficiency for possible muon decay electrons. Particle identification (PID) cuts 
are then applied to reject muon and 7r° events. Several improvements have been 
made to the neutrino oscillation data analysis since the initial data was published^, 
including an improved background estimate, an additional fiducial volume cut that 
greatly reduces the background from events produced outside the tank (dirt events), 
and an increase in the data sample from 5.579 x 10^° POT to 6.462 x 10^° POT. A 
total of 89,200 neutrino events pass the initial selection, while 1069 events pass the 
complete event selection of the final analysis with E^^ > 200 MeV, where E^^ is the 
reconstructed neutrino energy. 



3.4- Neutrino Oscillation Signal and Background Reactions 

Table[T]shows the expected number of candidate z/g CCQE background events with 
E^^ between 200 - 300 MeV, 300 - 475 MeV, and 475 - 1250 MeV after the complete 
event selection of the final analysis. The background estimate includes antineutrino 
events, representing < 2% of the total. The total expected backgrounds for the three 
energy regions are 186.8 ± 26.0 events, 228.3 ± 24.5 events, and 385.9 ± 35.7 events. 
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Table 1: The expected number of events in the 200 < < 300 MeV, 300 < E^^ < 475 MeV, and 
475 < < 1250 MeV energy ranges from all of the significant backgrounds after the complete 

event selection of the final analysis. Also shown are the expected number of ly^ CCQE signal events 
for two-neutrino oscillations at the LSND best-fit solution. 



Process 


ZUU — oUU 


oOO — 4/0 


/I 1 o txn 
4:10 — izoO 




y.u 


1 V /I 


ii. / 




6.1 


4.3 


6.4 


NC TfO 


103.5 


77.8 


71.2 


NC A iV7 


19.5 


47.5 


19.4 


Dirt Events 


11.5 


12.3 


11.5 


Other Events 


18.4 


7.3 


16.8 


z/g from /i Decay 


13.6 


44.5 


153.5 


Vf. from Decay 


3.6 


13.8 


81.9 


z/g from Decay 


1.6 


3.4 


13.5 


Total Background 


186.8 ±26.0 


228.3 ±24.5 


385.9 ±35.7 


LSND Best-Fit Solution 


7± 1 


37 ±4 


135 ± 12 



respectively. For — »• oscillations at the best-fit LSND solution of Am? = 1.2 
eV^ and sin^ 29 = 0.003, the expected number of Ue CCQE signal events for the three 
energy regions are 7 events, 37 events, and 135 events, respectively. 

3.5. Updated Neutrino Oscillation Results 

Fig. O shows the reconstructed neutrino energy distribution for candidate z/g data 
events (points with error bars) compared to the MC simulation (histogram) while 
Fig. [To] shows the event excess as a function of reconstructed neutrino energy. Good 
agreement between the data and the MC simulation is obtained for E^^ > 475 MeV; 
however, an unexplained excess of electron-like events is observed for E^^ < 475 
MeV. As shown in Fig. [TUl the magnitude of the excess is very similar to what is 
expected from neutrino oscillations based on the LSND signal. Although the shape of 
the excess is not c onsistent w ith simple two-neutrino oscillations, more complicated 
oscillation models I ^ l 5 | 6 | 7 | 8 |9t sterile neutrino decay ^ have shapes that may be 
consistent with the LSND signal. A test of the sterile neutrino decay model can be 
performed by searching for the decay D'^ —>■ fi'^i'h, where the heavy sterile neutrino 
i^h has a mass around 500 MeV. 

Table [2] shows the number of data, background, and excess events for different 
E^^ ranges, together with the excess significance. For the final analysis, an excess 
of 128.8 ± 20.4 ± 38.3 events is observed for 200 < E^^ < 475 MeV. For the entire 
200 < E^^ < 1250 MeV energy region, the excess is 151.0 ± 28.3 ± 50.7 events. As 
shown in Fig. [TTl the event excess occurs for Ey^g < 400 MeV, where E^ig is the 
visible energy. 
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Figure 9: The MiniBooNE reconstructed neutrino energy distribution for candidate Ve data events 
(points with error bars) compared to the Monte Carlo simulation (histogram). 
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Figure 10: The event excess as a function of E^^ . Also shown are the expectations from the best 
oscillation fit (sin^ 20 = 0.0017, Arn^ = 3.14 eV^) and from neutrino oscillation parameters in the 
LSND allowed region. The error bars include both statistical and systematic errors. 
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Tabic 2: The number of data, background, and excess events for different E^^ ranges, together with 
the significance of the excesses in neutrino mode. 
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Significance 
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Data 
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Background 
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Significance 
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Data 
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Significance 
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Figure 11: The event neutrino excess as a function of E^is for E^^ > 200 MeV. Also shown 
are the expectations from the best oscillation fit (sin^ 20 = 0.0017, Am^ = 3.14 eV^) and from 
neutrino oscillation parameters in the LSND allowed region. The error bars include both statistical 
and systematic errors. 



Figs. [12] and US] show the event excess as functions of and cos(6') for events 
in the 300 < E'^^ < 475 MeV range, where is determined from the energy and 
angle of the outgoing lepton and 6 is the angle between the beam direction and the 
reconstructed event direction. Also shown in the figures are the expected shapes 
from z/gC —>■ e'X and z/gC e^X charged-current (CC) scattering and from the NC 
tt" and A — > A^7 reactions, which are representative of photon events produced by 
NC scattering. The NC scattering assumes the energy spectrum, while the CC 
scattering assumes the transmutation of i/^ into Ve and Ve-, respectively. As shown in 
Table [3l the values from comparisons of the event excess to the expected shapes 
are acceptable for all of the processes. However, any of the backgrounds in Table [3] 
would have to be increased by > 5cr to explain the low-energy excess. 



3. 6. Initial Antineutrino Oscillation Results 



The same analysis that was used for the neutrino oscillation results is employed for 
the initial antineutrino oscillation results ^1^. Fig. UM shows the estimated neutrino 
fiuxes for neutrino mode and antineutrino mode, respectively. The fiuxes are fairly 
similar (the intrinsic electron-neutrino background is approximately 0.5% for both 
modes of running), although the wrong-sign contribution to the fiux in antineutrino 
mode (~ 18%) is much larger than in neutrino mode (~ 6%). The average plus z/g 
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Figure 12: The neutrino event excess as a function of for 300 < < 475 MeV. 
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Figure 13: The neutrino event excess as a function 0/003(6') for 300 < E^^ < 475 MeV. 
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Table 3: The values from comparisons of the neutrino event excess and cos{9) distributions 
for 300 < < 475 MeV to the expected shapes from various NC and CC reactions. Also shown is 
the factor increase necessary for the estimated background for each process to explain the low-energy 
excess and the corresponding number of sigma. 



Process 


X^{cose)/9 DF 


DF 


Factor Increase 


NC TfO 


13.46 


2.18 


2.0 (6.8a) 


A ^ iV7 


16.85 


4.46 


2.7 (18.4a) 




14.58 


8.72 


2.4 (15.3a) 


UeC e+X 


10.11 


2.44 


65.4 (41.0a) 



energies are 0.96 GeV in neutrino mode and 0.77 GeV in antineutrino mode, while the 
average z/^ plus z/^ energies are 0.79 GeV in neutrino mode and 0.66 GeV in antineu- 
trino mode. Also, as shown in Fig. [151 the estimated backgrounds in the two modes 
are very similar, especially at low energy. Fig. [16] shows the expected antineutrino 
oscillation sensitivity for the present data sample corresponding to 3.4E20 POT. The 
two sensitivity curves correspond to threshold neutrino energies of 200 MeV and 475 
MeV. 

The initial oscillation results for antineutrino mode are shown in Table [Hand Figs. 
[T7I through [T9l It is remarkable that no excess (— 0.5±7.8±8.7 events) is observed in 
the low-energy range 200 < E^^ < 475 MeV. In order to understand the implications 
that the antineutrino data have on the neutrino low-energy excess. Table [5] shows the 
expected excess of low-energy events in antineutrino mode under various hypotheses. 
These hypotheses include the following: 

• Same a: Same cross section for neutrinos and antineutrinos. 

• 7T^ Scaled: Scaled to number of neutral-current vr" events. 

• POT Scaled: Scaled to number of POT. 

• BKGD Scaled: Scaled to total background events. 

• CC Scaled: Scaled to number of charged- current events. 

• Kaon Scaled: Scaled to number of low-energy kaon events. 

• Neutrino Scaled: Scaled to number of neutrino events. 

Also shown in the Table is the probability (from a two-parameter fit to the data) 
that each hypothesis explains the observed number of low-energy neutrino and an- 
tineutrino events, assuming only statistical errors, correlated systematic errors, and 
uncorrelated systematic errors. A proper treatment of the systematic errors is in 
progress; however, it is clear from the Table that the "Neutrino Scaled" hypothe- 
sis fits best and that the "Same a", "POT Scaled", and "Kaon Scaled" hypotheses 
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Figure 14: The estimated neutrino fluxes for neutrino mode (top plot) and antineutrino mode (bottom plot). 
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Figure 15: The estimated backgrounds for the neutrino oscillation search in neutrino mode (top plot) and antineutrino 
mode (bottom plot). The tt", A — > Nj, intrinsic Ve/i^e, external event, and other backgrounds correspond to the green, 
pink, light blue, blue, and yellow colors, respectively. 
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Figure 16: The expected antineutrino oscillation sensitivity at 90% CL for the present data sample 
corresponding to 3.4E20 POT. The two sensitivity curves correspond to threshold energies of 200 
MeV (red curve) and ^15 MeV (black curve). 
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Table 4: The number of antineutrino data, background, and excess events for different ranges, 
together with the significance of the excesses in antineutrino mode. 



Event Sample 


Final Analysis 


200 - 475 MeV 




Data 


61 


Background 


61.5 ±7.8 ±8.7 


Excess 


-0.5 ±7.8 ±8.7 


Significance 


-0.04a 


475 - 1250 MeV 




Data 


61 


Background 


57.8 ±7.6 ±6.5 


Excess 


3.2 ±7.6 ±6.5 


Significance 


0.3cr 


475 - 3000 MeV 




Data 


83 


Background 


77.4 ±8.8 ±9.6 


Excess 


5.6 ±8.8 ±9.6 


Significance 


OAa 



are strongly disfavored. It will be very important to understand this unexpected 
difference between neutrino and antineutrino properties. 

The antineutrino data were also fit for oscillations in the energy range 475 < 
E^^ < 3000 MeV, assuming antineutrino oscillations but no neutrino oscillations. 
The antineutrino oscillation allowed region is shown in Fig. [201 At present, the 
oscillation limit is worse than the sensitivity. The best oscillation fit corresponds to 
Am^ = 4.4 eV^, sin^ 26* = 0.0047, and a fitted excess of 18.6 ± 13.2 events, which is 
consistent with the LSND best-fit point of Am^ = 1.2 eV^ = 0.003, and an 

expected excess of 14.7 events. With the present antineutrino statistics, the data are 
consistent with both the LSND best-fit point and the null point, although the LSND 
best-fit point has a better {x^ = 17.63/15 DF, probability = 30%) than the null 
point (x^ = 22.19/15 DF, probability = 10%). 

3. 7. MmiBooNE NuMI Results 

Neutrino events are also observed in MiniBooNE from the NuMI beam^^. The 
NuMI beam, as shown in Fig. [211 differs from the Booster neutrino beam (BNB) in 
several respects. First, the NuMI beam is off axis by 110 mrad, whereas the BNB 
is on axis. Second, neutrinos from NuMI travel ~ 700 m, compared to ~ 500 m for 
neutrinos from the BNB. Also, the NuMI beam has a 6% contribution from electron- 
neutrinos and a 14% contribution from antineutrinos, while the BNB percentages 
are 0.5% and 2%, respectively. Fig. [221 shows the estimated neutrino flux at the 
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Figure 17: T/ie MiniBooNE reconstructed antineutrino energy distribution for candidate Ve data 
events (points with error bars) compared to the Monte Carlo simulation (histogram). 



Table 5: The expected excess of low-energy events in antineutrino mode under various hypotheses 
for 3.4E2O POT. Also shown in the Table is the probability (from a two-parameter fit to the data) 
that each hypothesis explains the observed number of low-energy neutrino and antineutrino events, 
assuming only statistical errors, correlated systematic errors, and uncorrelated systematic errors. 



Hypothesis 


^ of u Events 


Stat. Err. 


Cor. Syst. Err. 


Uncor. Syst. Err. 


Same a 


37.2 


0.1% 


0.1% 


6.7% 


7r° Scaled 


19.4 


3.6% 


6.4% 


21.5% 


POT Scaled 


67.5 


0.0% 


0.0% 


1.8% 


BKGD Scaled 


20.9 


2.7% 


4.7% 


19.2% 


CC Scaled 


20.4 


2.9% 


5.2% 


19.9% 


Kaon Scaled 


39.7 


0.1% 


0.1% 


5.9% 


Neutrino Scaled 


6.7 


38.4% 


51.4% 


58.0% 
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Figure 18: The MiniBooNE reconstructed antineutrino energy distribution for candidate 9^ data 
events (top) and the excess number of events (bottom) as a function of reconstructed neutrino en- 
ergy for the present antineutrino data sample corresponding to 3.4E20 POT. Also shown are the 
expectations from the best oscillation fit and from oscillation parameters in the LSND allowed re- 
gion. 
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Figure 19: The excess number of candidate events (data minus Monte Carlo expectation) as a 
function of visible energy for the present antineutrino data sample corresponding to 3.4E20 POT. 
Also shown are the expectations from the best oscillation fit and from oscillation parameters in the 
LSND allowed region. 



MiniBooNE detector from the NuMI beam, while Fig. [23] compares the neutrino 
fluxes from the BNB and NuMI beams. 

Figs. [21] and [25] show the comparison between data events (points with error 
bars) and the MC simulation (histogram) for i/^ CCQE candidate events and Ve 
CCQE candidate events, respectively. Although the systematic errors are presently 
large, the data are observed to be systematically low for z/^ CCQE candidate events 
and systematically high for z/g CCQE candidate events. Updated results should be 
available soon with three times the data sample and with reduced systematic errors 
by constraining the normalization to the sample. 

The NuMI data analysis is currently directed toward examining the low-energy 
region and searching for neutrino oscillations. This will complement the analysis 
done with MiniBooNE using neutrino and anti-neutrino BNB data, but with different 
systematic errors. It is worth noting that the NuMI z/g CCQE sample has a very 
different composition when compared to the BNB neutrino CCQE sample. The 
BNB Vf. CCQE sample originates mostly from decays of pions and muons and contains 
a large fraction of z/^-induced mis-identified events. On other hand, the NuMI z/g 
CCQE sample is produced mostly from the decay of kaons and contains a dominant 
fraction of intrinsic z/g events. The analysis will be done by forming a correlation 
between the z/^ CCQE and z/g CCQE samples and by tuning the prediction to the 
data simultaneously. The result is that common systematics cancel, and this might 
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Figure 20: The antineutrino oscillation allowed region in the energy range 475 < E'^^ < 3000 MeV 
for the present antineutrino data sample corresponding to 3.4E20 POT. Also shown are the best 
oscillation fit (^w? — AA eV^ , sin^ 29 = 0.0047, corresponding to an excess of 18.6 ± 13.2 events) 
and the LSND best-fit point (Am^ = 1.2 eV^ , sin^ 20 = 0.003, corresponding to an excess of 14.7 
events). 
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Figure 21: The NuMI beam. 
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Figure 22: The estimated neutrino flux at the MiniBooNE detector from the NuMI beam. 
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Figure 23: A comparison between the BNB and NuMI neutrino fluxes. 
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Figure 24: The comparison between data events (points with error bars) and the MC simulation 
(histogram) for NuMI-induced CCQE candidate events. 
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Figure 25: The comparison between data events (points with error bars) and the MC simulation 
(histogram) for NuMI-induced CCQE candidate events. 



reveal something important about the nature of the z/g sample. 



3.8. MiniBooNE Disappearance Results 

MiniBooNE has also searched for z/^ and disappearance Fig. [26] shows 
the sensitivities and limits at 90% CL for z/^ disappearance (top plot) and z/^ dis- 
appearance (bottom plot). The stars on the plots show the best fits in each case: 
Am^ = 17.50 eV^, sin^ 26^ = 0.16, and = 12.72/14 DF for i/^ disappearance and 
Am^ = 31.30 eV^, sin^ 29 = 0.96, and = 5.43/14 DF for z/^ disappearance. The 

values for no disappearance oscillations are = 17.78/16 DF and = 10.29/16 
DF, respectively. Improved disappearance sensitivities are expected with the joint 
SciBooNE/MiniBooNE analysis, which should be completed soon. Note, however, 
that the joint SciBooNE/MiniBooNE analysis will not be nearly as powerful as the 
joint BooNE/MiniBooNE analysis due to the small size of SciBooNE and the larger 
systematic errors from SciBooNE's different detector technology. 

4. BooNE 

The BooNE experiment involves building a second detector at a cost of ~ $8M 
along the BNB at FNAL at a closer distance of ~ 200 m. With two detectors, 
many of the systematic errors will cancel, as the neutrino flux varies as 1/r^ to 
good approximation, so that a ratio of events in the two detectors will provide a 
sensitive search for z/g and z/g appearance and z/^ and z/^ disappearance. Furthermore, 
by comparing the rates for a NC reaction, such as NC scattering or NC elastic 
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Figure 26: The sensitivities and limits at 90% CL for disappearance (top plot) and disappear- 
ance (bottom plot). 
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scattering, a direct search for sterile neutrinos can be made. An even cheaper option 
would be to move the MiniBooNE detector to a different location at a cost of only ~ 
$4M. If the MiniBooNE detector were moved to a distance of 200 m from the neutrino 
source, then the event rate would increase by a factor of ~ 7 due to the dependence of 
the neutrino flux on distance. An additional advantage of moving MiniBooNE is that 
MicroBooNE could then move into the original MiniBooNE enclosure and, therefore, 
save the expense of building a new MicroBooNE enclosure. In either case, after less 
than a year of running, the comparison of the event rates at the two locations will 
determine whether the low-energy excess observed by MiniBooNE was due to neutrino 
oscillations. In addition, i/^ and z/^ disappearance will be searched for with high 
sensitivity in the Am^ > 0.1 eV^ mass region, and LSND antineutrino oscillations can 
be tested directly by searching for z/g appearance. By comparing neutrino oscillations 
to antineutrino oscillations, BooNE will be able to search for CP and CPT violation 
in the lepton sector at short baseline (Am^ > 0.1 eV^). For the sensitivities discussed 
below, it is assumed that the near detector will run for ~ 1E20 POT in both neutrino 
mode and antineutrino mode. 

4.1. Fluxes and Event Rates 

This section gives a detailed comparison of the expected neutrino fluxes at the 
near (200 meter) and far (541 meters) positions. In the Booster neutrino beam 
(BNB), the primary beam is produced by the 8 GeV Fermilab's rapid-cycling (15Hz) 
booster accelerator, which produces 1.6 fis batches of protons each containing around 
4.5 X 10^^ protons. 

At that primary proton energy, there are only four significant species of neutrinos: 
z/^ and (~ 99.5%), and a small contamination (~ 0.5%) of u^. and z/g. There are 
two primary parent components to the fluxes: neutrinos from charged pion decays 
and neutrinos from kaon decays. The component dominates the z/^ spectrum 
above neutrino energies of 2.5 GeV, where a clear break is observed in the slope of 
the energy spectrum. The z/^ spectra are mainly from charged pion decay, and the z/g 
and z/g spectra are composed of two parts, muon decays and kaon decays. 

The standard MiniBooNE Geant4 based beam simulation and decay program 
packages were used to generate fluxes^'''*. Those packages include the transport of 
muon polarization (neglecting g — 2 precession effects) and appropriate form factors 
in leptonic kaon decays. The primary production of pions by 8 GeV protons was 



measured by the HARP experiment!^^ and is used as input in the simulation, while 
secondary interactions in the beam line are handled by standard Geant4 physics 
packages. 

The fluxes shown here represent the spectrum of neutrinos that intersect a sphere 
of radius 610.6 cm, positioned at either the near or far location. The fluxes are 
"unoscillated" and therefore have only z/^(z/^)and z/e(z/e) components. No matter 
effects in propagating the neutrinos to the detector are included, as they are expected 
to be small in the standard, 3-generation, active neutrino model (Sz/M). 
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Figs. [271 EHl EH and [30] show the fluxes for the four neutrino species at the near 
and far locations, for both neutrino mode and antineutrino mode. Table [6] gives the 
same fluxes, integrated over neutrino energy, while Table [7] gives the average neutrino 
energy in each case. 

In neutrino mode, the u^^ flux near/far ratio is 7.5. Most of the near far ratios are 
between 7.0 and 8.0. Another characteristic of the near/far flux comparisions is that 
the average energy of the neutrinos in the near position is between 5 and 10 percent 
less than the corresponding average energy in the far position. This lower energy is 
expected since the near detector has a larger angular acceptance with respect to the 
neutrino target. 

Fig. [3T] shows the energy distribution for CCQE events at the near (1.0 x 10^° 
POT) and far (6.462 x 10^° POT) locations for neutrino mode. The spectral differences 
are again due to the larger angular acceptance of the near detector. That larger decay 
angle translates to lower neutrino energies in the near detector, typically ~ 10% lower 
in the 200/541 meter comparison. This extrapolation is relatively straight forward as 
the angular divergence of the daughter neutrinos in the decays is much larger than 
the angular divergence of the decaying mesons themselves. For example, even at 3 
GeV, daughter neutrinos from pion and kaon decays will have opening angles of ~ 
50 mrad and ~ 150 mrad, respectively, while the allowed angular divergence of the 
beam tunnel is only ~ 20 mrad. 

Because of the nearly complete overlap in decay particle phase space that con- 
tributes to neutrinos in the near and far positions, we expect that uncertainties in the 
flux prediction will largely cancel when comparing the two event rates from the near 
and far positions. As systematic errors introduced by uncertainties in the detector 
efficiency and neutrino cross section will also largely cancel, the comparison of the 
two positions will allow a much-needed, accurate measurement of non-Si^M neutrino 
oscillation effects in the Am^ range of 0.1-10 eV"^. 

Table 6: Integrated fluxes per POT for the various species of neutrinos at the near and far positions, 
for both neutrino mode and antineutrino mode. 



Fluxes z//(cm2 POT) 





V mode 


V mode 


v species 


Near 


Far 


Near 


Far 




7.49 X 10-8 


1.03 X 10' 


-8 


8.12 X 10-9 


1.08 X 10- 


-9 




5.20 X 10-9 


6.52 X 10- 


10 


4.30 X 10-*^ 


5.77 X 10' 


-9 




4.50 X 10-1° 


5.74 X 10- 


11 


9.5 X 10-" 


1.34 X 10- 


11 




4.61 X 10^11 


6.00 X 10- 


12 


2.00 X 10-1° 


2.53 X 10- 
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Jf-.2. Possible Scenarios for a Near Detector 

The MiniBooNE detector has operated at a location of 541 meters from the 
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Table 7: Average neutrino energies for the various species of neutrinos at the near and far positions, 
for both neutrino mode and antineutrino mode. 



Average v energies (MeV) 
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Figure 27: The and u^^ fluxes at both the near and far locations in neutrino mode. 



Booster Neutrino Beam (BNB) target since September 1, 2002. The primary pur- 
pose of the experiment was to search for the transmutation, or oscillation, of muon 
neutrinos into electron neutrinos as they travel the ~ 525 meters to the detector. 
The BNB was designed to produce a nearly pure beam of z/^, which provides an ideal 
setting to look for excess Ve events. While the most sensitive neutrino oscillation ex- 
periments are two detector systems, which afford a comparison of similar detectors at 
two different distances, MiniBooNE was built as a single detector system in order to 
reduce costs. It was felt that the systematic error incurred by not building a second 
detector could be overcome by using internal measurements in the single detector. 
As were not expected to oscillate significantly, it was planned to use their rate as 
a normalization for v^. interactions, thus constraining backgrounds to u^. events from 
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Figure 28: The and fluxes at both the near and far locations in antineutrino mode. 
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Figure 29: The Ve o,nd i>e fluxes at both the near and far locations in neutrino mode. 



oscillations. 

The MiniBooNE proposal foresaw that a second detector at a different distance 
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Figure 30: The and fluxes at both the near and far locations in antineutrino mode. 




ure 31: The distribution of energies for reconstructed CCQE events, at the near location with 
X 10^° POT and at the far location with the current 6.462 x 10^'' POT. 
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would be required to ascertain the nature of the signal, if a significant signal were ob- 
served in the single-detector setup. A second detector, located at a different distance 
from the BNB target could potentially remove the large systematic errors that would 
complicate the interpretation of the MiniBooNE data. 

There are several possible routes to a more precise measurement which need to be 
considered: a new detector could be constructed at a near position; similarly, a new 
detector could be constructed at a far position; and, it was recently realized that the 
MiniBooNE detector could be relocated to a new position at either a near or a far 
location. Each of those possibilities has advantages and disadvantages. 

4.2.1. 

Near or Far? 

The choice between constructing a near detector at ~ 200 meters and a far detector 
at ~ 1000 meters can be made based on expediency. For Am^ < 2 eV^, a near detector 
will not see a large signal directly but can be used to accurately measure the expected 
backgrounds to any possible oscillation signal in the far detector. Those backgrounds, 
in both appearance and disappearance measurements, can be measured at ~ 7-8 times 
the rate that the MiniBooNE detector accumulated data. Thus a sample of neutrinos 
with statistics equivalent to MiniBooNE's existing data set of 6.462 x 10^° POT will 
only require ~ 1.0 x 10^° POT and yield an ~ 5cr result. An identical far detector, 
on the other hand, would also yield an ~ 5cr result with ~ 1.0 x 10^° POT, however 
the signal would only be ~ 20 events on top of a background of ~ 16 events. An 
unsatisfying, ambiguous result could occur with such low statistics. 

4.2.2. 

Moving MiniBooNE 

Relocating the MiniBooNE detector to a near position ~ 200 meters from the 
target, shown in Fig. [32l will likely be the least costly option. There are a number of 
potential advantages to this approach. The detector is already built so that the cost 
of constructing a new detector is avoided. The neutrino flux at 200 meters will be ~ 
7-8 times larger, so that the time needed to accumulate a data sample equivalent to 
the existing sample will be less than one year. 

The comparison of measurements taken by the same detector operated at two 
locations will be free of systematic errors associated with neutrino-nucleus cross sec- 
tions and detector response. The comparison should clearly reveal the nature of the 
MiniBooNE excess. The 541/200 meter data comparison will also allow MiniBooNE 
to measure z/^ and disappearance at the few percent level. 

The detector relocation could take place in two ways: transporting the existing 
detector and electronics to the new location, or building a new detector at the new 
location using parts from the existing detector. In both cases the MiniBooNE dec- 
tector would be drained of oil and the oil stored in rail cars or a separate oil storage 
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Figure 32: An arial view of Fermilab showing a possible location of a near detector at ~ 200 meters. 
The blue arrow indicates the direction of the booster neutrino beam. 

tank. With the mineral oil removed, the detector weighs only ~ 80 tons, and could 
be transported whole to the near site. Mobile cranes with 80 ton Ifting capacities are 
readily available commercially at a reasonable cost. 

4.2.3. 

A New Detector 

An alternative to relocating the existing MiniBooNE detector and its electronics 
is to construct a new sperical tank at the near position and remove PMTs by hand 
from the old tank. The new tank would be constructed by repeating the constuction 
effort made for MiniBooNE. In either case, one would re-use the existing, albeit old, 
LSND electronics currently used by MiniBooNE. The cost of either case is estimated 
to be ~ $4M. 

The most desirable, and most costly, option is to construct an entirely new detector 
at the near location. This would require more time because new electronics would have 
to be developed, a new oil mineral supplier found, and new phototubes purchased. 
The lead times on those items would be about two years. The cost for that effort is 
estimated to be ~ twice that of moving the existing detector, ~ $8M. 

It is not yet understood how the systematic error in detector response will translate 
between the old MiniBooNE detector and the newly constructed detector, since it 
will have different oil, PMTs, and electronics. Nevertheless, choosing to construct an 
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Table 8: The expected excess of low- energy events (200 < < 475 MeV) in antineutrino mode 
(1E21 POT) for various hypotheses, assuming a 2.5% systematic error and assuming that the neu- 
trino low-energy excess is correct. Also shown is the significance of the excesses. 



Hypothesis 


Expected Excess of v Events 


Significance 


Same a 


111.6 ±13.5 ±4.6 


7.8a 


tt'' Scaled 


58.2 ±13.5 ±4.6 


4.1a 


POT Scaled 


202.6 ± 13.5 ±4.6 


14.2a 


BKGD Scaled 


62.8 ±13.5 ±4.6 


4.4a 


CC Scaled 


61.2 ±13.5 ±4.6 


4.3a 


Kaon Scaled 


119.2 ± 13.5 ±4.6 


8.4a 


Neutrino Scaled 


20.2 ±13.5 ±4.6 


1.4a 



entirely new detector would allow for simultaneous operation of both the near and 
far detector and eliminate any fear, unfounded or not, that the neutrino beam had 
changed in character. 

Ideally, the near BooNE detector would have the same dimensions as the Mini- 
BooNE detector in order to reduce systematic uncertainties. However, another possi- 
bility would be to build a smaller detector (~ 8 m diameter) at a lower cost (~ $4M) 
if systematic errors were estimated to be sufficiently small. 

4-3. Testing the Low-Energy Excess 

BooNE will be able to determine whether the low-energy excess is due to neutrino 
oscillations and will be able to test various hypotheses by comparing the low-energy 
excesses in neutrino and antineutrino modes. If the low-energy excess is due to 
background, then the near detector will observe the same relative excess as the far 
detector. If the excess is due to neutrino oscillations at low Am^, then no low-energy 
excess will be observed in the near detector and the current low-energy excess in 
the far detector, assuming a 2.5% systematic error, will equal to 128.8 ± 20.4 ± 10.4 
events (5.6a). For testing various hypotheses. Table [S] shows the expected excess 
of low-energy events (200 < E^, < 475 MeV) in antineutrino mode for 1E21 POT, 
assuming a 2.5% systematic error and assuming that the neutrino excess is correct. 
Also shown is the significance of the excesses. By comparing the measured excess in 
the antineutrino data to the expected excesses, the different hypotheses can be shown 
to be either consistent with data or ruled out. 

4.4- and Ue Appearance 

The sensitivities for z/g and appearance were obtained by using the MiniBooNE 
Monte Carlo simulation, assuming statistical errors with the expected MiniBooNE 
statistics (6.5E20 POT in neutrino mode and 1E21 POT in antineutrino mode) and a 
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full error matrix with correlated and uncorrelated systematic errors. Also, we assume 
2E20 POT in the near detector, equally divided between neutrino and antineutrino 
modes. Fig. [33] shows the estimated sensitivity for Vf, appearance for Ey > 200 
MeV. Although simple two-neutrino oscillations have already been ruled out as an 
explanation of the LSND signal, it is interesting that the full LSND region is covered 
at the approximately 5a level. Therefore, we would be able to determine whether 
or not the MiniBooNE low-energy excess is due to a more complicated oscillation 
mechanism at the ~ 1 eV^ scale. Fig. [M] shows the estimated sensitivity for z/g 
appearance, where we assume that the error matrix is the same as for neutrinos. The 
sensitivity is worse than the z/g appearance sensitivity due to the lower statistics and 
higher wrong-sign background in antineutrino mode; however, BooNE will still be 
able to cover the full LSND region at 90% CL and provide a direct test of LSND 
antineutrino oscillations. 

4-5. Ufj, and Disappearance 

The sensitivities for and disappearance were obtained by using the Mini- 
BooNE Monte Carlo simulation, assuming statistical errors with the expected Mini- 
BooNE statistics (6.5E20 POT in neutrino mode and 1E21 POT in antineutrino 
mode) and a full error matrix with correlated and uncorrelated systematic errors. 
Also, we assume 2E20 POT in the near detector, equally divided between neutrino 
and antineutrino modes. Fig. [351 shows the estimated sensitivity for z/^ disappearance 
for Ey > 200 MeV. A sensitivity of ~ 3% at 90% CL is obtained for Am^ ~ 1 eV^. In 
order to see how a signal would appear. Fig. [361 shows the allowed region for z/^ disap- 
pearance at the global antineutrino best-fit point from reference Amf^ = 0.915 
eV^ and sin^ 2^^^ = 0.35. Figs. [37] and [3H] show the correponding limits and al- 
lowed regions for z/^ disappearance, assuming no z/^ disappearance and the same error 
matrix as for neutrinos. The sensitivity is slightly worse than the z/^ sensitivity 
due to the lower antineutrino statistics and the ~ 1/3 wrong-sign z/^ component in 
antineutrino mode. A difference between z/^ and z/ ^ disappearance would be evidence 
for CPT violation or effective CPT violation I^SEl 

4-6. Sterile Neutrino Search 

If z/^ or Un disappearance is observed, then the NC tt" and NC Elastic reactions 
can be used to determine whether the disappearance is due to oscillations into active 
or sterile neutrinos. Oscillations into sterile neutrinos will result in a suppression of 
events in the far detector, while oscillations into active neutrinos will result in no 
suppression. Due to the high statistics of the NC vr^ and NC Elastic event samples, 
the statistical error will be small compared to the systematic error. The sensitivity 
at 90% CL for oscillations into sterile neutrinos at Am^ ~ 1 eV^ is estimated to be 
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Figure 33: The estimated BooNE sensitivity for Ve appearance. 
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Figure 34: The estimated BooNE sensitivity for Ve appearance. 
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Figure 35: The estimated BooNE sensitivity for disappearance. 
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Figure 36: The BooNE allowed region for v^^ disappearance at the global antineutrino best-fit point 
Am?, = 0.915 eV^ and sin^ 26l.„ = 0.35. 
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Figure 37: The estimated BooNE sensitivity for disappearance. 
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Figure 38: The BooNE allowed region for v^^ disappearance at the global antineutrino best-fit point 
Am?, = 0.915 eV^ and sin^ 26l.„ = 0.35. 
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sin^ 26^^^ ~ 3% for neutrinos and sin^ 26^^ ~ 5% for antineutrinos. 

5. Other Experiments 

The MINOS experiment has measured neutrino osciUations at the atmospheric 
scale with z/^ disappearance and has begun to look at disappearance 133. It is 
interesting to note that the MINOS antineutrino data so far are consistent with the 
antineutrino 3+1 model of reference and consistent with i?^ disappearance at the 
LSND scale. Fig. [39] shows the MINOS z/^ event rate as a function of reconstructed 
neutrino energy compared to the Monte Carlo expectation. Although the statistics 
are low, the MINOS event rate is suppressed above 10 GeV, where there should 
be almost no suppression due to atmospheric neutrino oscillations. Fig. HQ] shows the 
MINOS z/^ disappearance allowed region, which is consistent with oscillations at the 
~ 1 eV^ scale. 

The only other experiment or proposal that is capable of addressing these physics 
objectives at the ~ 1 eV^ scale is a Letter of Intent to refurbish the CERN PS 
neutrino beam and build two liquid argon detectors . However, the proposed 
BooNE experiment, with the existing BNB, should be able to obtain results prior to 
the CERN experiment. 

6. Cost Estimate 

Table [H] shows a breakdown of the cost estimate for constructing a second BooNE 
detector in a near location. The estimate is based on the MiniBooNE construction 
costs. The total estimated cost is $7.3M, including contingency (~ 30%) and escala- 
tion (3% per year). The BooNE construction is assumed to start in 2010 and last for 
3 years. The estimated cost for moving MiniBooNE to a near location is ~ $4M. An 
additional advantage of moving MiniBooNE is that the MicroBooNE detector could 
then be moved into the original MiniBooNE enclosure, thereby saving the expense of 
building a new MicroBooNE enclosure. 
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Figure 39: The MINOS event rate as a function of reconstructed neutrino energy compared to 
the Monte Carlo expectation. Although the statistics are low, the MINOS event rate is suppressed 
above 10 GeV, where there should be almost no suppression due to atmospheric neutrino oscillations. 



Item 


Cost ($K) 


Tank and support structure 


1065 


PMT's 


1759 


Electronics /DAQ 


512 


Oil 


1429 


Calibrations 


412 


Miscellaneous 


198 


Engineering & Construction 


1894 


Total 


7269 



Table 9: A breakdown of the cost estimate for constructing a second BooNE detector in a near 
location, including contingency and escalation. 
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